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An exact integro-differential equation for the conformal map from the unit circle to the boundary 
of an evolving cavity in a stressed 2-dimensional solid is derived. This equation provides an accurate 
description of the dynamics of precursors to fracture when surface diffusion is important. The 
solution predicts the creation of sharp grooves that eventually lead to material failure via rapid 
fracture. Solutions of the new equation are demonstrated for the dynamics of an elliptical cavity 
and the stability of a circular cavity under biaxial stress, including the effects of surface stress. 

PACS numbers: 46.50.+a, 62.20.Mk,81.40.Np 
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The process of rapid fracture of solids which have failed 
to sustain stress is a poorly understood subject Q ; elas- 
ticity theory does not suffice to describe this process, 
since plastic deformations occur at the most interest- 
ing "process zone" where the actual fracture is taking 
place 0. The displacement field is not the only relevant 
field, and even if the stress field is given everywhere, it 
is not known how the fracture propagates: there exist 
complex interactions with sound waves and maybe other 
fields 0>Q • On the whole it is not obvious how to achieve 
a self consistent theory. 

The situation is much clearer when one studies slow 
processes that may precede rapid fracture. In particular 
we will discuss in this Letter "precursors" to fracture. 
Namely, the dynamics of stress driven deformations of 
cavities (or free surfaces) in solids 0, IE S S S EH EH 
Such deformations are expected to lead, even- 
tually, to the creation of deep grooves which then suffer 
such large stresses that the solid fails via rapid fracture. 
The aim of this Letter is to derive and demonstrate a 
new equation of motion for the conformal map from the 
unit circle to the evolving stress-driven deformed cavity 
in 2-dimcnsional solids. This equation offers an accu- 
rate description of the slow evolution of the precursors 
to failure, until the moment that rapid fracture can be 
sustained. In contrast to all previous treatments we in- 
clude both surface energy and surface stress and show 
that our equation is well posed. 

To set up the problem imagine a 2-dimensional elas- 
tic medium with a hole of an arbitrary shape, whose 
boundary z(s) is parametrized by the arc-length variable 
s. Boundary conditions at infinity load the medium, in- 
ducing a displacement field u(r) with a strain tensor ejk, 
related to the stress by 
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Here v is the Poisson ratio and E the Young modulus. 
Under stress there begins a process of surface diffusion 
which deforms the boundary, with dynamics determined 



by the velocity v n normal to the boundary, 
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where D is a diffusion coefficient that depends on the 
material and the temperature, and [i is the chemical po- 
tential at the boundary [l6|. The chemical potential is 
obtained from the change in total energy associated with 
an infinitesimal variation of the interface [|| 0] : 
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where 7, k and (3 are a the surface energy, the mean cur- 
vature and surface stress respectively. S stands for the 
deformation energy S = j tij&ij /2 The notation t 
and n stands for "tangent" and "normal" to the interfce, 
defined precisely in Eqs.©. Mo and C are the reference 
chemical potential and a material parameter that we can 
take as zero and unity respectively. Note that we differ 
from H, licj in taking into account the surface stress. 
Ref. [I3 incorporated all this physics but performed sta- 
bility analysis of the fiat interface only. 

To evolve the cavity we need to compute then the 
chemical potential on its boundary, consistent with the 
evolving stress field in the medium. We will assert that on 
the slow time scale of surface diffusion the elastic medium 
is in equilibrium, i.e. £\ dk <Jjk = for all j. The 
general solution of these equations in 2-dimensional is 
given by 
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where the so called Airy potential x fulfills the bihar- 
monic equation A 2 % = 0. The general solution of this 
equation is written in complex notation, with z = x + iy, 

z = x — iy, as 



X(z,z) = 3? z(f>(z) +i>(z) 



(5) 



where 4>{z) and ip(z) are any pair of analytic functions, 
to be determined from the boundary conditions. 



2 



Consider the boundary of the cavity. With a the an- 
gle between the tangent and the s-axis at z(s), define 
derivatives in the tangent and normal directions: 



d t — cos(a) d x + sin (a) d y 
d n = — cos(a) d y + sin(a) d x 



(6) 



The surface stress now must be balanced by the normal 
component of the stress 0, 0] : 
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On the other hand the mixed derivatives vanish since 
there is no restoring force along the boundary, 



dtnX = otn = 0nt = on the crack 



(8) 



Using these boundary conditions we can evaluate the 
chemical potential /i on the boundary [In principle /i 
has terms that cannot be written in terms of Tra alone, 
however these terms add to zero as a consequence of 
dn&nn = —dt<J n t and da n t / ds = on the boundary] 



1 



E 



(Tra) 2 



(3d n Tra — finTraj — 7K 



(9) 



Using the fact that Ad 2 x/ dzdz — a xx + a yy — Ott + <Jnn 
we can immediately read from Eq. (JSJ), 

Tr<r = 4St[4>\z)] ■ (10) 

Thus to compute [i and its derivatives and advance the 
cavity we only need to determine the function 4>{z). The 
boundary conditions Q and (jSJl are expressed in terms 
of <f){z) and tjj(z) by using Eqs. © We note that 

(cos a + % sin a)(pt — id n ) = 29z and on the boundary 

1 (3k 
dtd z x = -(cosa + isma)d tt X = ~^^t z > ( n ) 

where we identify dtz(s) as the unit vector tangent 
to the boundary. Writing the mean curvature as 
K = dfz(s)/idtz(s), this condition reads dt{dzX + 
i(3dtz(s)/2} = 0. Thus the boundary condition on the 
interface [n| is 



4>(z(s))+z(s) <f>'(z(s))+ip(z(s)) = -i(3d t z(s) + K , (12) 

where ip(z) = ^'(z) and if is a constant (that can be 
chosen zero with impunity). 

For the boundary conditions at infinity we consider 
biaxial loading: 
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It was shown in 19] that if the integral over the boundary 
of the RHS of Eq. H12|l is zero, then the finiteness of 
the stresses at infinity implies that the analytic functions 
have a Laurent expansion of the form 



<f>(z) = (f>iz + Ei=o^-* 2 

1p(z) = IpiZ + Ysi^^-i*' 



(14) 




FIG. 1: Evolution of an ellipse under a biaxial load. Shown 
is the change in the radius Ap t (6) = |$(e, t)| — |$(e, 0)| (in 
dimensionless units) relative to its initial value. The plot is 
given at dimensionless times 10 4 x t = 4, 30, 40, 58.9, 71.8.77.0 
and 82.3. In the insert we show the analogous evolution for a 
cavity without surface stress. 



The freedoms that we have in determining the Airy func- 
tion JHJ) allow us to choose 4>q — and cj>i real. Then 
using the boundary conditions (|13|l at infinity, we find 
4>i = <T o/2, "01 = 0. To proceed consider now a con- 
formal map &(uj,t) which maps the exterior of the unit 
circle e = exp(i#) to the exterior of boundary z(s). Our 
central aim in this Letter is to derive an equation of mo- 
tion of this map. The normal velocity v n is related to 
the conformal map as follows: denote by h — n x + in y 
and z{s) = x{s) + iy(s) the unit vector normal to the 
boundary, and the position of the boundary respectively 
. Then 



v n (s) = x(s)n x + y(s)n y = $t 



dz{s) ■ 
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dt 



(15) 

(Note that below we use e and exp(i8) interchangeably). 
The tangent (t) and normal unit vector are given in terms 
of the conformal map by 



gg(e) _ d9 gg(e) 
ds 89 
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Using this in Ea. ll5|) and remembering that d&(e l6 )/dt 
d t § + <&'e lB i9 t we derive: 



$'(e)| 



(17) 
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We rewrite this equation, 



d>'(e) 



= u„ + iC 



(18) 



with an unknown imaginary part C. Multiplying the last 
equation by e<j>'/|<j>'| we get the equivalent equation 



<9f$ = e$'(e) 
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and unknown C". This equation, valid on the interface, 
can be analytically continued outside the unit circle. We 
need to choose C such that the term in the parenthesis 
is an analytic function, removing the freedom in C . This 
also fixes the parameterization which was so far arbitrary. 
All this is achieved with the Poisson integral formula, 
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The equation, being analytic, must have analytic solu- 
tions which provide the dynamics of the conformal map. 
In practice, the Poisson integral formula is best expressed 
in terms of the Fourier components of the function. Given 
a real function on the unit circle 
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the analytic function outside the circle whose real part 
on the unit circle is g is 



G = 



271 



d6< 
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T g(e ie ) = a + 2 Y a n o 



(22) 



Thus finally the equation of the conformal map reads 



9 t $ = uj&(w)G(u>) . 



(23) 



The conformal map itself (which is univalent) has a Lau- 
rent expansion of the form 



*(w,t) = Fi(t) u + F (t) + 



n=l 
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By taking derivatives with respect to t and lo and sub- 
stituting back in Ea. (|23[) we reduce the dynamics to an 
infinite set of ordinary differential equations for the Lau- 
rent coefficients Fi(t). All that remains is to compute 
the function 17(e) that is given in terms of v n (e) and the 
conformal map. To compute the normal velocity u n (e) 
first compute Trtr on the boundary. Rewrite Ea. i|12|) in 
terms of the conformal map at time t: 



0'($(e,i))+^($(e,t)) 
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FIG. 2: The normal velocity as a function of the angle 6 on 
the unit circle. The times are the same as in Fig. The 
insert shows the analogous figure for a cavity without surface 
stress. 



According to Eq.(j2U the functions in Eq. 1)25(1 read 

0($(e,t)) - fcFi(t)e + 0o(e) 

lK$(e,t)) = ^Fi(t)e + ^(e) , (26) 

where the only unknowns now are the expansion coeffi- 
cients of the functions 4>o and tpo, 
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At this point we need to compute the two Fourier trans- 
forms 



= E »i 



n— — 00 



|$'(e,t)| 



00 

E ^ 

71= — CO 



(28) 



substitute them in to Eq. |25jl . and obtain a system 
of equations for the coefficients u_„ and In fact 

upon the substitution one finds that the determination 
of the function 4>(u) is independent of ip(uj), requiring 
only negative powers of e. 

Having found the function (f>(u>) we employ Eq. (|1()[1 to 
compute the trace of the stress tensor anywhere on the 
boundary. The curvature K is: 



1 $" -n 
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Taking the second derivative of the chemical potential 
we have v n (e). With this function at hand we can re- 
turn to Eq. (|21|l and step the equations for the Laurent 
coefficients, finding the new conformal map, etc. 
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an initial elliptical hole and (ii) the stability analysis of 
a circular hole. For (i) our initial conformal map is 



$(w,t = 0) = Fi(0)u + 



F-i(0) 



(30) 
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FIG. 3: The curvature at the tip as a function of time for an 
ellipse under compression and extension, and for zero surface 
stress. Dotted line: the curvature predicted by the linear 
stability analysis. 



To demonstrate the efficiency of this procedure, and 
the interesting predictions it offers for the dynamics of 
cavities in stressed solids, we present (i) the evolution of 



where we chose Fi(0) = 1, F_i(0) = 0.01. For material 
parameters we took (1 — v 2 )jE = 1, 7 = 0.2, D = 1 
and (3 — 0.1, all in dimensionless units. We study the 
evolution of the initial ellipse under biaxial load with 
(To = 1. In Fig. ^ we show the time evolution of p t = 
|$(e, t)| — |$(e, 0)|. In the inset we show the analogous 
dynamics when (3 = 0, 7 = 0.1. We see that the tip 
of the ellipse is advancing at the expense of two dips 
that develop on its side, but this effect is more dramatic 
for (3 = as seen in the inset. In Fig. 2 we show the 
normal velocity, with the inset for the same parameters 
as in Fig.l. Finally, in Fig. [21 we show the curvature at 
the tip. When the stress at the tip reached the Griffith 
criterion the material would yield by rapid fracture |20j |. 
For (ii) we start with 

$(«,*) = Rio + Y,fn{t)u;- n ; \f n \ « 1 Vn , (31) 

ri=l 

which maps the unit circle (e = e %e ) to a wavy shaped 
circle of radius R. Using our equation l|23l for the con- 
formal map, and Eq. (|25(l which determines the stress, 
linearizing in small /„, it is a straightforward calculation 
to obtain the stability eigenvalues in the form 
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Contrary to the calculation in [Tl| our problem is well 
posed, the cubic term (n 3 ) is negative. The analytic re- 
sult l|32|) is novel, showing precisely what load is needed 
to destabilze a circular cavity. 

In summary, we stress that surface stress term in fi 
may matter: it introduces for sharp tips a competitive 
term proportional to r~ 3 l 2 in contrast to the terms pro- 
portional to 1/r of the stress energy and the curvature, 
surface stress removes the degeneracy of compression and 
extension, as we see in Eq. (|32H . In the nonlinear regime 



the grooves start to form faster for compression. These 
and related issues will be elaborated further in a future 
publication. 
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